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Abstract 

It is argued considering multiple production in the deep inelastic scattering 
{Sj . kinematics that the very high multiplicity events are extremely sensitive to the 

$h ! low-x partons density. 

^ T J * 

\q : 1 Introduction 

The role of soft color partons in the high energy hadron interactions is the mostly 
t— i \ intriguing modern problem of particles physics. So, the collective phenomena and sym- 

metry breaking in the non-Abelian gauge theories, confinement of colored charges and 
the infrared divergences of the pQCD are the phenomena just of the soft color particles 
domain. 

It seems natural that the very high multiplicity (VHM) hadron interaction, where 
the energy of created particles is small, should be sensitive to the soft color particles 
Ph! densities. Indeed, the aim of this paper is to show that even in the hard by definition 

deep inelastic scattering (DIS) the soft color particles role becomes important in the 
VHM region. 

The standard (mostly popular) hadron theory considers perturbative QCD (pQCD) 
at small distances (in the scale of A ~ 0.2 Gev) as the exact theory. This statement is 
confirmed by a number of experiments, e.g. the DIS data, the QCD jets observation. 
But one should have in mind that the pQCD has finite range of validity since the non- 
perturbative effects should be taken into account at distances larger then 1/A. 

It is natural to assume building the complete theory that at large distances the non- 
perturbative effects exceed the perturbative onesQ. In result the pQCD contributions 
becomes screened by the non-perturbative interactions. 

But exist another, more natural, possibility. It should be noted here that the pQCD 
running coupling constant ce s (q 2 ) = l/61n(g 2 /A 2 ) becomes infinite at q 2 = A 2 and we do 
not know what happens with pQCD if q 2 < A 2 . There is a suspicion [[| that at q 2 ~ A 2 
the properties of theory changed so drastically (being defined on new vacuum) that even 
the notions of pQCD is disappeared. This means that pQCD should be truncated from 
below on the 'fundamental' scale A. It seems evident that this infrared cut-off should 
influence on the soft hadrons emission in the VHM domain. 

So, it is important try to rise predictability of pQCD in the 'forbidden' area of 
large distances. For this purpose one should split experimentally the perturbative and 



1 The corresponding formalism was described e.g. in Q. 
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non-perturbative effects at the large distances to check out quantitatively the role of 
soft color partons. One must realize for this purpose highly unusual condition that the 
non-perturbative effects must be negligible even if the distance among color charges is 
high. 

The non-perturbative effects lead to the strong polarization of QCD vacuum and, in 
result, to the color charge confinement. This vacuum is unstable @ against creation of 
real quark- ant iquark (qq) pairs if the distance among charges became large. This pairs 
are captured into the colorless hadrons and just emission of this 'vacuum' hadrons is the 
mostly important non-perturbative effect. Therefore, generally, the number of hadrons 
n, produced even in the hard DIS process, n ^ n p , where n p is the number of qq pairs 
created 'perturbatively'. 

Notice, if the kinetic energy of colored partons is small, i.e. is comparable with 
hadron masses, creation of 'vacuum' hadrons should be negligible. Just this is the VHM 
process kinematics: because of the energy-momentum conversation law, produced (final- 
state) partons can not have high relative momentum and, if they was created at small 
distances, the production of 'vacuum' hadrons will be negligible (or did not play im- 
portant role). Therefore, if the 'vacuum' channel is negligible, the pQCD contributions 
only should be considered [Q. 

The aim of this article is to show that n — > oo unavoidably leads to 'low-x' domain. 



2 DIS kinematics in the VHM domain 

To describe the hadron production in pQCD terms the parton-hadron duality is assumed. 
This means that the multiplicity, momentum etc. distributions of hadron and colored 
partons are the same. This reduce the problem practically on the level of QED. 

Let us consider now n particles (gluons) creation the DIS. We would like to calculate 
D ab (x,q 2 ;n), where 

J2Dab(x,q 2 ;n) = D ab (x,q 2 ). (2.1) 

n 

As usual, let D ab (x, q 2 ) be the probability to find parton b with virtuality q 2 < in the 
parton a of ~ A virtuality, A >> A and a s (X) << 1. We always may chose q 2 and x so 
that the leading logarithm approximation (LLA) will be acceptable. One should assume 
also that (1/x) >> 1 to have the phase space, into which the particles are produced, 
sufficiently large. 

Then D ab (x, q 2 ) is described by ladder diagrams. From qualitative point of view this 
means approximation of Markovian process of random walk over coordinate ln(l/x) and 
time is In In \q 2 \. LLA means that the 'mobility' ~ ln(l/x)/ lnln \q 2 \ should be large 

ln(l/x) » lnln q 2 /X 2 . (2.2) 

But, on other hand |J, 

ln(l/x) «ln|g 2 /A 2 |. (2.3) 

The leading, able to compensate smallness of a s (X) << 1, contributions give inte- 
gration over wide range A 2 << k 2 << —q 2 , where k 2 > is the 'mass' of real, i.e. 
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time-like, gluon. If the time needed to capture the parton into the hadron is ~ (1/A) 
then the gluon should decay if k 2 » A 2 . This leads to creation of (mini)jets. The mean 
multiplicity fij in the QCD jets is high if the gluon 'mass' [ A: j is high: Inn,,- ~ ^J\n(k 2 / A 2 ). 

Rising multiplicity may (i) rise number of (mini)jets v and/or (ii) rise the mean value 
mass of (mini)jets < \ki\ >. We will see that the mechanism (ii) would be favorable. 
This is consequence of the Markovian character of considered process. 

But rising mean value of gluon masses \ki\ decrease the range of integrability over 
ki, i.e. violate the condition ( |2.2| ) for fixed x. One can remain the LLA taking x — > 0. 
But this may contradict to ( |2.3|) , i.e. in any case the LLA becomes invalid in the VHM 
domain and the next to leading order corrections should be taken into account. 

Noting that the LLA gives main contribution, that the rising multiplicity leads to 
the infrared domain, where the soft gluons creation becomes dominant. 



3 Preliminary notes 

First of all, neglecting the vacuum effects, we introduce definite uncertainty to the 
formalism. It is reasonable to define the level of strictness of our computations. Let us 
introduce for this purpose the generating function T ab (x, q 2 \ z): 

D ab (x, q 2 ; n) = ^— <j> -^T ab {x, q 2 ; z). (3.1) 

Having large n the integral may be calculated by the saddle point method. The smallest 
solution z c of the equation 

d 

n = z— lnT ab (x, q 2 ;z) (3.2) 
oz 

defines the asymptotic over n behavior: 

D ab {x, q 2 ; n) oc e -****>Mri . (3.3) 
Using the statistical interpretation of z c as the fugacity it is natural to write: 

\nz c {x,q 2 -n) = Cab{X ' q \^ . (3.4) 
n a b{x,q 2 ) 



Notice that the solution of eq.(p.2|) z c (x,q 2 ;n) should be the increasing function of n. 
At first glance this follows from positivity of all D ab (x, q 2 ;n). But actually this assumes 
that T ab (x,q 2 ; z) is a regular function of z at z = 1. This is a natural assumption 
considering just the pQCD predictions. 
Therefore, 

D ab (x, q 2 - n) oc e~*^& Cah{x '*''' n) . (3.5) 

This form of D ab (x, q 2 ; n) is useful since usually C ab (x, q 2 ; n) is the slowly varying func- 
tion of n. So, for Poisson distribution C ab (x,q 2 ;n) ~ Inn. For KNO scaling we have 
C ab (x, q 2 ;n) = const, over n. 
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We would like to note, that neglecting effects of vacuum polarization we introduce 
into the exponent so high uncertainty assuming n ~ n p that it is reasonable perform 
the calculations with the exponential accuracy. So, we would calculate 

- fx ab (x, q 2 ; n) = In D °*( x ' q L^j = _ U C ab (x, q 2 ; n)(l + 0(1/ n)) (3.6) 

The n dependence of C ab (x,q 2 ;n) defines the asymptotic behavior of fi ab (x,q 2 ;n) 
and calculation of its explicit form would be our aim. 



4 Correlation functions 

Considering particles creation in the DIS processes one should distinguish correlation 
of particles in the (mini)jets and the correlations between (mini)jets. We will describe 
the jet correlations. One should introduce the v jets creation inclusive cross section 
<f>P"(ki, k 2 , k v ; q 2 , x), where fcj, i = 1,2,..., n are the jets 4-momentum. Having <3>j, 
we can find the correlations function N^(ki, k 2 , k v \ g 2 , x), where (r) = r±, r v and 
r i — (?> 9) defines the sort of created color particle. It is useful introduce the generating 
functional 

F ab (q 2 , x;w)=J2 I dV n {k) f[ w r *(k t ) fc 2 , k n , q 2 , x)\* , (4.1) 

n J i=l 

where a ab is the amplitude, dQ n (k) is the phase space volume and w Ti (ki) are the arbi- 
trary functions. It is evident, 

F ab (q 2 ,x;w) \ w=1 = D ab (q 2 , x) . (4.2) 

The inclusive cross sections 

k 2 , k u ; g 2 , x) = n F ab (q 2 , x; w) \ w=1 . (4.3) 
The correlation function 

NP(k u h, K- g 2 , x) = f[ ^ F ab (q 2 , x; w) \ w=1 . (4.4) 

We will use this definitions to find the partial structure functions D ab (q 2 ,x;n). 
It will be useful to introduce the Laplace transform over variable ln(l/x): 

F ab (q 2 ,x;w)= [ jL(lyf*tfj ;w ) (4.5) 
JRcj<o 2m x 

The expansion parameter of our problem a s ln(— g 2 /A 2 ) ~ 1. By this reason one 
should take into account all possible cuttings of the ladder diagrams. So, calculating 
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D ab (q 2 , x)in the LLA all possible cuttings of sceleton ladder diagrams is defined by the 
factor §: 

- {Vf G r Vf\ , (4.6) 

i.e. the cutting line may get not only through the exact Green function G r {kf) but 
through the exact vertex functions T^ b (qi, qi+i, hi) also (qf , q 2 +1 are negative). We have 
in the LLA that 

\ 2 , ^ 2 ^ ^ 2 ^ ^ 2 

A << -q\ « -q l+1 « -q 

and 

X < Xi + \ < Xi < 1. 

Following to our approximation, see previous section, we would not distinguish the way 
as cut line go through the Born amplitude 



{(Tf) 2 G r } 



We will simply associate w r lma.f°) to each rung of the ladder. 

Considering the asymptotics over n the time-like partons virtuality ki ~ —qf/yi 
should be maximal. Here yi is the fraction of the longitudinal momentum of the jet. 
Then, slightly limiting the jets phase space, 

In kl = In |g m | 2 (1 + 0(ln(l/x)/|g m | 2 )). (4.7) 

In result, introducing useful in the LLA variable Tj = ln(g, 2 /A 2 ), where ot s (q 2 ) = 
1//3t, (3 = (11A/3) — (2n//3), we can find following set of equations: 

d 

r-K-fab(q 2 ,j;w) = J2 ( P r acU) wr ( T )fab(q 2 ,x;w), (4.8) 

<J ' c,r 

where 

flcij) = Vac(j) = / —X>Pac(x) (4.9) 

Jo X 

and P ac (x) is the regular kernel of the Bethe-Salpeter equation 0. At w — 1 this 
equation is the ordinary one for D ab (q 2 ,x). 

We will search the correlation functions from eq. (|4.8 ) in terms of Laplace transform 
n^ v (ki, k2, k v \ q 2 ,j). Let us write: 

(4.10) 

Inserting ( |4.10| ) into ( |4.8| ) and expanding over (w — 1) we find the sequence of coupled 
equation. 

Omitting the cumbersome calculations, we write in the LLA that 

t 2 , r u ; q 2 ,j) = d aci {j, ri)^ C2 (j)4 2C3 (j, t 2 ) ■ ■ ■ (pZc„ +1 (j) d e v+1 b(j, t„+i). 

(4.11) 
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One should take into account the conservation laws: 

ri • r 2 • • ■ r u+1 = t, n <t 2 < ... < T u+ i < r. (4.12) 

Computing the Laplace transform of this expression we find (t 1; r 2 , r v ; q 2 , x). 
The kernel G? b(j, r) was introduced in ( 4.11| ). Let us write it in the form: 

d ab {j,T) = Y.°^-^\ (4.13) 

(T=± T 



where 
and 



= "a ~ <Pgg, dqg = v a - (p n , d a qg = <p gq , <F gq = <p qg (4.14) 

V *=\ \fn + ¥99 + ° [(<Pqq ~ <P BB Y ~ ^ f^ qg <f gq ] ^ } . (4.15) 

If x « 1, then (j — 1) << 1 are essential. In this case ||, 

This means the gluon jets dominance and 

n 9 gg = <p 99 + 0(l). (4.17) 

One can find following estimation of the two-jet correlation function: 

<r(ri,r 2 ; J) r) = 0(max{(r 1 /r)^ ! (r 2 /r)^ ) (r 1 /r 2 )^}}. (4.18) 

This correlation function is small since in the LLA T\ < t 2 < r. This means that the 
jet correlations becomes high iff the mass of correlated jets are comparable. But this 
condition shrinks the range of integration over r and by this reason one may neglect the 
'short-range' correlations among jets. Therefore, as follows from (|4.10| ), 

fab{q 2 ,j;w) = d gg {T,j)exp\ip gg [ ^-w 9 (r')l (4.19) 



/t r 

We will use this expression to find the multiplicity distribution in the DIS domain. 



5 Generating function in the DIS kinematics 

To describe particles production one should replace iy r Ima^ b ) on ^Imaf), where w r n is 
the probability of n particles production, 

= (5-1) 

n 

Having v jets one should take into account the conservation condition n\ +n 2 + ... +n u = 
n. By this reason the generating functions formalism is useful. In result one can find 
that if we take (PD 

w 9 = w 9 (r,z), w 9 (r,z) \ z=1 = 1, (5.2) 
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then fab(<l 2 ,j',w) denned by (|4.19|) is the generating functional of the multiplicity dis- 



tribution in the l j representation'. In this expression w 9 (t } z) is the generating function 
of the multiplicity distribution in the jet of mass \k\ = Xe T ^ 2 . 
In result, see fl4.5|) , 

F ab (q 2 ,x;w) oc / -^(1/ z )V» w(t '* ) (5.3) 
JRej<0 2m 

where 

FT df 1 

u(t,z)= / —w 9 (r',z). (5.4) 

Jtq t 



Noting the normalization condition (J57 

lu(t, z — 1) = lnr. (5.5) 



The integral ( |5.3| ) may be calculated by steepest descent method. It is not hard to 
see that 

j ~ j c = 1 + {4A^(r, z)/ ln(l/a;)} 1/2 (5.6) 

is essential. Notice that j — 1 << 1 should be essential we find, instead of the constraint 
Q, that 

w(t,z) « ln(l/x). (5.7) 

In the frame of this constraint, 



F ab (q 2 , x; w) oc exp U^Nu(t, z) ln(l/x) j . (5. 



Generally speaking, exist such values of z that j c — 1 ~ 1. This is possible if u;(r, z) 
is a regular function of z at z = 1. Then z c should be the increasing function of n and 
consequently cu(t, z c ) would be the increasing function of n. Therefore, one may expect 
that in the VHM domain j c — 1 ~ 1. 

Then j ~ 1 + uj(t, z)j ln(l/x) would be essential in the integral (|5.3|) . This leads to 
following estimation: 

F ab (q 2 ,x;w) oce-"^. 

But this is impossible since F ab (q 2 ,x;w) should be the increasing function of z. This 
shows that the estimation ( |5T8| ) has a finite range of validity. 

Solution of this problem with unitarity is evident. One should take into account 
correlations among jets considering the expansion ( f4.1Up . Indeed, smallness of n^" may 



be compensated by large values of Hi w ri (T i} z) in the VHM domain. 

6 Conclusion 

We can conclude that our LLA is applicable in the VHM domain till 

w{t, z) « ln(l/x) «r = \n(-q 2 /\). (6.1) 
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The mean multiplicity of gluons created in the DIS kinematics 



d 



n g {r,x) = — In F ab {q 2 ,x;w) | z=1 = uj x {t)^ AN\n{l/x)/ lnr » wi(r), (6.2) 
where 

Wl (r)= f—n^r) (6.3) 
and the mean gluon multiplicity in the jet nj(r) has following estimation ||: 

lnnj(r) ~ v^" (6-4) 



Inserting ( |Q| ) into ( |6T3| 
Therefore, noting ( |2.3|) , 



umt) = n,- r /Vr. 



(r,x) ~ nj(T)yj4Nhi(l/x)/ThiT « %(r). (6.5) 



This means that the considered 't-channel' ladder is important in the narrow domain of 
multiplicities 

n ^ fig « fly (6.6) 

So, in the VHM domain n » n g one should consider 

(i) The ladder diagrams with small number of rungs; 

(ii) To take into account the malty-jet correlations assuming that increasing multiplicity 
leads to the increasing number of rungs in the ladder diagram. 

To choose one of this possibilities one should consider the structure of oo(t, z) much more 
carefully. This will be done in the consequent paper. 

We can conclude that the VHM domain multiplicities production unavoidably de- 
stroy the ladder LLA. To conserve this leading approximation one should choose x — > 
and, in result, to get to the multi-ladder diagrams, since in this case a s ln(— q 2 /X 2 ) ~ 1 
and a s ln(l/x) ~ 1. Such theory was considered in |7[. 
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